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Abstract 

Using heavy quark symmetry we construct an energy sum rule relating the 
form factor of a heavy flavored meson £(w) to the light quark energy distri- 
bution. We find that the current available data for £(w) is consistent with a 
broad energy distribution rather than with that of a pure quark/anti-quark 

bound state, indicating a large nonvalence content. 
12.38Lg, 13.20Jf, 14.40Jz 
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There is growing interest in the physics of heavy quark systems stimulated by the recog- 
nition [T|J^] of the important role of the heavy quark symmetry of the QCD lagrangian in the 
infinite quark mass limit. Interest has also been fostered by recent measurements of many 
heavy hadron decay modes [[| which provide more stringent challenges of our understanding 
of the interplay between strong and weak interactions [|J. In the heavy quark mass limit, 
m — > oo, the spectrum of heavy mesons containing one heavy and one light valence quarks, 
become degenerated and independent of the polarization state of the light degrees of free- 
dom 0. Furthermore, an effective decoupling of the dynamics of a heavy quark reduces the 
problem to that of a cloud of light quarks and gluons in an external colored source |§[5]||. 
These two facts imply definite relations between heavy meson matrix elements of various 
currents that involve heavy quark field operators. In particular it has been shown that in 
the limit of an exact heavy quark symmetry all form factors associated with such matrix 
elements are determined by a single, universal function £(w), called the Isgur-Wise form 
factor The function £(«;) is usually introduced through a matrix element of a vector 

heavy quark current evaluated between two J p = 0~ heavy meson states, 

(p'|¥(o)7"*(o)b) = U(W + pY + f-(W ~pY> (!) 

with p' 2 = M' 2 , p 2 = M 2 and t = (p' — p) 2 . In the m, M, M' — > oo limit it is useful to 
introduce the velocity transfer, t 2 /m 2 = (V — v) 2 . Defining w = v ■ v' — 1 — t/2m 2 the form 
factors f±(t), which are not independent reduce to a single function of w, 

M' + M 

f+(t) -> ; = £M->eM, 

M - M' 

There have been many attempts to calculate the universal function £(w) ari d a ^ so 

the 0(l/m) corrections to f± ||T0| , [TT|] . As already pointed out, due to the decoupling of 
the heavy quark, is primarily determined by the dynamics of light degrees of freedom 
thereby considerably simplifying the problem. This point is implicit in the constituent quark 
model calculations || in which £(«;) is determined by a noncovariant wave function of a light 



spectator quark. Another approach is based on a QCD sum rule analysis of the vacuum 
three-point correlator ||. The standard interpolation method is used to relate the spectral 
representation of the correlator to the contribution from the ground state, thus to In 
such an approach the decoupling of the heavy quark makes £(w) depend on the propagator 
of the light quark in an external field generated by the heavy quark. 

In this paper we are interested in studying the structure of the distribution of light 
degrees of freedom in the heavy meson, and propose an approach in which we relate £ (w) to 
the nondiagonal correlator of two currents, i.e. with only one of the meson states in Eq. ([!]) 
being interpolated. The correlator is then evaluated between the vacuum and the other 
heavy meson state. It is in a sense an intermediate approach between the constituent quark 
model and a three-point function QCD sum rule calculation, and it will become clear that 
such an approach will allow us to relate, in a covariant way, the universal function to 
gauge invariant matrix elements describing light quark and gluon content of a heavy meson. 
The general form of a correlator which we shall be dealing with in the rest of this paper is 
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1.1 



T(p,q) = tj dz^MTOfyoA—m, (3) 



where \p,i) is a heavy meson state with momentum p, (p 2 = M 2 ) containing a light 
valence quark with flavor i, and the operators 0,Oi are given by 0(z) = ^/(z)T^/(z), 
Oi(z) = ^/(z)T'tlJi(z). Here \& and ipi are the heavy and light quark field operators and 
T and T' are some combination of Dirac gamma matrices, respectively. We start with 
the phenomenological analysis of the correlator T(p,q). The 2-dimensional plane in vari- 
ables q\ 2 i (<?i,2 = p/2 ± q), is conventionally parameterized in terms of u, A and q 2 , i.e. 
q\ 2 = (l±uj)q 2 ±\ + ^. For fixed u and A, with \u\ < 1 and uX > p 2 /A, T(p,q) = T(u,X,q 2 ) 
satisfies a dispersion relation in the variable q 2 [12| with 

5 l3 ImT(X,u;,q 2 ) = \j dze* [<p, *|O(|)P,-(-|)|0> + (p, i|0,(-|)0(f )|0) 

Inserting a complete set of intermediate states with masses M n in Eq. gives, 

SijT(X,u,q 2 ) 



(4) 



I ]_ — |— TYVtTTt j 

S M 2 _ A _ p2/ /4_ (1 + ^ )g2 _ i€ (P^|O(0)l W .gln)Kgln|O 3 -(0)|0) 9 ; w =( P /2+ 9 „) a =A^ 



+A~P 2 /4 _ 



M 2 + A _ p2/4 _ (1 _ w)g2 *l°i(°)i n > Q2n)(n, g 2n |O(0)|0) gL=(p/2 _ gn)2=M 2. 



(5) 



where s TO j n = min (—^7^, J 4 ) . With our choice of the quark content for the operators 
O, Oj the first term in the RHS of Eq. (|5|) has contributions from resonances containing 
one heavy and one light valence quark, while the second term involves states with two 
heavy quarks and does not contribute to £(w). Therefore, we shall choose the variables 
u,X,q 2 in such a way that the second term is nonleading in the 1/m expansion. This 
corresponds to the choice A > p 2 /4, (u > 0) and q 2 ~ p 2 /2(l + lj) which leads to a 
suppression by a factor of 0(l/m) of the second term with respect to the first term in 
Eq. (0). Furthermore it can be shown that the (1 — u) dependence is also nonleading so 
that in the leading order analysis we set A > p 2 /A and u = 1. It is also more convenient 
to use, in the infinite heavy quark mass limit |13| , [T^| , the spectral representation in the 
energy variables instead of the invariant mass. With the resonance fixed at q\ n = M 2 the 
momentum transfer is t — 2m 2 (l — w) = (p — qi n ) 2 = —A + p 2 /4 and the correlator will be 
evaluated for q\ = (m + E q ) 2 ~ m 2 + 2mE q . Denoting the energy of an excited state by E n , 
M 2 ~ m 2 + 2mE n , the correlator in the variables E q and w for u = 1 and A > p 2 /4 is given 
by 

5 ijT ( W , E q ) = E 2m(E n -E q -ie) " (6) 

In the following we restrict our analysis to two choices of the operators O, Oj given by two 
sets of the T, V matrices; Set 1 : T — 7^, V = 75, Set 2 : T = 1, r" = 75. We will explicitly 
show that to leading order in the 1/m expansion the two sets lead to a unique sum rule for 
the form factor £(«;). For \E q — E n \ >> a standard assumption is that the contributions to 
the sum in Eq. (^|) from the ground state and from higher resonances can be parameterized 
by a single state with energy E and zero width together with a continuum starting at 



with threshold energy E c For two T sets the leading contribution to the phenomenological 
correlator as m — > oo is given by 



&Al:T*{w,E q ) = ±(lp + q)i> 



MM 



Set 2 : T(w, E q ) 



% 

—m 

2 



(E q -E + ie 
f h (l + w)^w) 



+ T c (i, E q , E c ) 



T c (t, E q , E c 



(7) 



_(E q -E + ie) 

where f\ is a ground state heavy meson decay constant, E c is the effective continuum 
threshold energy, and T c is the contribution from the continuum whose explicit form will 
be discussed later. In deriving the second equation above we have used (p'|\I/(0)\I/(0) |p) = 
m(l + w)£(w), which follows from Eq. ([[]) taking the derivative of the vector current in the 
M' — M limit. 

We shall now discuss the theoretical description of the correlators. For E q >> Eq ~ 
Aqcd the heavy quark at the intermediate state in the correlator in Eq. ([|) is off its energy 
shell roughly by E q and the singularity of the perturbative heavy quark propagator dominates 
the spacetime behavior of the operator product in Eq. (H). In a fixed gauge, the leading 



perturbative contribution is given by [18 



(8) 



where S(z) is the perturbative heavy quark propagator and the trace is taken over the spinor 
indices. The nonperturbative part of the propagator leads to a 0(l/m) correction and is 
not taken into account in the leading order analysis. For the two sets of operators discussed 
earlier the correlators are given by 



T»(p,q) = if h Jd A z J 



d 4 k 



-i(q—k)z 



(2vr) 4 k 2 - m 2 + it 
+ k fi m(j)_(z 2 , z • p) + k u m(pTufi(z 2 , z ■ p) 

e 



+ p^)m(j)(z 2 , z ■ p) 



T(p, q) = if h d z 



d 4 k P -i(<}-k)z 



(2vr) 4 k 2 - m 2 + ie 



(k ■ p + m )4>{z ,z-p)+m 4>-{z ,z-p) 



(9) 



where 



4>\z 2 ,z-p) 



t(z 2 ,z-p) = — 

Jh 



2 p z 2 

1 ,z 



m 



2p 2 



4>Tvu{z 2 , Z ■ P) 



T^(0|^(|)^M75*(-|)|p,i), 



(10) 



and the normalization of the leading amplitude, <j)(z 2 , z ■ p), is given by <ft(z 2 , z ■ p) z =o — 1- 
In the infinite heavy quark mass limit the operators that depend on a spin structure of the 
light quark (i.e, \E r 7 M 75^j and ^75^) are degenerate and thus lead to the same matrix 
element. For this reason <p_ does not contribute in the leading order. Similarly it is seen 
that <pT is suppressed by a power of 1/m relative to <fi so only terms proportional to <fi will be 
retained in the leading order analysis. The matrix element which defines can be rewritten 



as 



! :-p) = ^e^(0|^(z)475*(0)b,*> = &'%(£,z L ) 



(11) 



with the subscript L indicating the magnitude of a 4-vector projection parallel to p M and 
subscript T indicating the remaining three components of a 4-vector perpendicular to p^. 
The structure of the light quark wave function, <f) q is revealed by expanding in a Taylor series 
at z — 0, 



\{z Tl z L , 



E 



{izl^/p 2 ) 



nl 



:)(o)+«+tOT+ 



(12) 



The moments (m) (0) can be related to local matrix elements with n longitudinal, D^, 
and 2m transverse, Dt, covariant derivative insertions, 



d» + igA^(z) (Tj. If we 
imagine that the heavy meson is composed of a single light quark bound by a chromoelectric 
potential of a static heavy quark, i.e. with no dynamical gluons, then Di = 8l + igA^zx) 
being z^ independent implies that (4>g)^ m '(0) with fixed m determine the energy (longitudinal 
component of the 4-vector) distribution of a single light quark. The heavy meson having 
a definite energy, Eq, requires the light quark to be on the energy shell as well, and as a 
consequence (0p(O) = (E /^ s /p I ) n or 

q (z 2 ,z L )=e^ E ^ q (z 2 ). (13) 
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The transverse amplitude <p q {z^) is the bound state wave function which has a nontrivial 
behavior for \zr\ ~ l/^o even in the static case and it approaches a plane wave solution 
only in a free case i.e. when A(zt) - ► 0. However the zl dependent amplitude is the same 
whenever <p q describes a bound state or a free particle. Since in general the vector potential 
A = A(zl, zt) is nonstatic one expects a sizable contributions from the gluon and/or qq sea. 
The z L dependence of A L implies that the spectral representation of <f> q (zi,) is a smeared 
distribution around E = Eq rather then being proportional to a delta function. Furthermore 
since the gluon field A^ = A^{z) should be in general p M independent, the division between 
the longitudinal and transverse wave functions is purely a matter of convenience. Also 
relevant is the scale governing the energy and momentum smearing intervals of the two 
wave functions which should be of the same order. Thus, despite the fact that in leading 
order the correlators in Eq. (^j are not sensitive to the transverse wave function, the above 
argument still permits conclusions on the size of the gluon momentum distribution from the 
analysis of the longitudinal wave function alone. 

Transforming to the momentum representation the correlators are given by 

% 3 . f dE <f) q {E) 



T"(w,E q )= l -f h (-p + qrJ o 



o 2-7T E q — Ew + ie ' 



T(w, E q ) = L f hm {l + w) f ^i E) (14) 

v ' q> T v ; Jo 2vr E q - Ew + ie K J 

with (j)q(E) = J dz L e iEzL 4>(z^ = 0,z L ) and J ^-<p q (E) = 1. The leading order theoretical 
correlators of Eq. ( |T4D are also used to define a continuum contribution T c introduced in 
Eq. (0). The standard assumption is that the spectral density of T c (t, E q , E c ) representing 
the sum over higher resonances can be replaced by the theoretical one at E > E c . From 



Eqs. (ffl) and ( 141) we obtain a single sum rule for £(w). 



dE <p q {E) £(w) + r dE <p q {E) ^ 



lo 2tt E q - Ew + ie E q - E + ie JEz 2ir E q - Ew + ie 
In order to reduce the contributions from the unknown higher order terms in the expansion 
in Eq. ( |T2"D and suppress the sensitivity of the sum rule to the phenomenological parameter- 
ization of the spectral density, a standard Borel transformation (E q —>■ T) is performed for 



Eq. (|15D leading to 

dE 



e~ wE °/ T + I d ±e~ wE l T \UE) - 27r5(E - E )] 
Jo Aix 



Z(w)e- E °/ T + f d 4e- wE/T UE). (16) 



As explained earlier in a simple quantum mechanical description in which a heavy meson 
contains a single quark orbiting around a static chromoelectric source, 4> q {E) = 2tt5(E — E ) 
and since E c > Eq, the sum rule automatically leads to a correct normalization, £(1) = 1, 
for all values of T. At first look the above result may seem dependent on the choice of 
the parameterization of the phenomenological spectral density. This is however not the 
case. If a more detailed parameterization were used, involving explicitly higher resonances 
instead of a smooth continuum approximation, due to orthogonality of the mass eigenstates 
such resonances would not contribute to the RHS of Eq. (|T6"D at w = 1 leaving £(1) = 1 
coming only from the ground state contribution. In a realistic situation 4> q (E) is a smeared 
distribution around E = E . Furthermore it follows from the analysis of the normalization, 
fH-*Li of the heavy-to-light meson matrix elements that the /h^l ~ 0(m~ 3 / 2 ) behavior [15 



requires 4> q (E) to vanish only linearly with E for E/Eq « 1 |fT6| . In order to be consistent 
with the phenomenological parameterization of the spectral density we shall assume that 
<p q {E) vanishes for E > E c . The value of the threshold energy has been obtained in Ref. jUTl 



and it varies in the range 2E < E c < 3E . We then use the following parameterization 
(fiq oc E{E C — E)exp(— (E — E ) 2 /E^) and study the form factor for different values 
of the width parameter Ew- For a given value of E c /Eq and Ew/Eq the Borel parameter 
(T/E ) is fixed by the normalization £(1) = 1. In Fig. 1 we plot our predictions for £(u>) 
for E c /E = 2.1,2.5,3 given by the set of dashed, solid and dotted curves respectively and 
for E w /E = .1 and 10 corresponding to upper and lower curves in each set (dashed, solid, 
dotted) respectively. As E w decreases the predictions are sensitive to both the detailed shape 
of the energy distribution of the light quark and to the value of the continuum threshold 
energy. On the other hand, for increasing Ey/ our predictions become weakly dependent 
on E c /Eq and Ew/Eq, and are bound from below by 1/w 2 behavior which in turn follows 
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from the linear behavior of <j) q (E) at small E. The data corresponding to the B — > D form 
factor (finite m) and extracted from the B — > Dlv\ with A = (tb/1.18) * (|V^,|/0.05) = 



1.11 have been taken from Ref. ||19|| . The comparison of our results with the experimental 
data shows that the £(w) form factor is inconsistent with the usually assumed peaking (8- 
type) approximation to the distribution amplitude corresponding to small Ew/E ) ^(J and 
suggests a rather broad, E w > E , energy distribution. This in turn implies large gluon 
and possibly sea quark amplitudes in a heavy meson and we conclude that a large fraction 
of the energy-momentum of the light degrees of freedom is distributed over the nonvalence 
components. Although our results rely on a comparison with heavy, but finite, quark mass 
data, both the 0(l/m) and perturbative, 0(a(E)) corrections, as shown in Ref. are 
expected to be very small not exceeding a few percent. 
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FIGURES 

FIG. 1. The Isgur-Wise form factor £(w). The dash-doted line is the 1/w 2 asymptotic behavior. 
The variuos theoretical curves are explainded in the text. 
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